DO STRANGE STARS EXIST IN THE UNIVERSE? 
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Abstract. Definitely, an affirmative answer to this question would have 
implications of fundamental importance for astrophysics (a new class of 
■ compact stars), and for the physics of strong interactions (deconfined 

phase of quark matter, and strange matter hypothesis). In the present 
work, we use observational data for the newly discovered millisecond X- 
ray pulsar SAX J1808.4-3658 and for the atoll source 4U 1728-34 to con- 
strain the radius of the underlying compact stars. Comparing the mass- 
radius relation of these two compact stars with theoretical models for 
both neutron stars and strange stars, we argue that a strange star model 
^ ! is more consistent with SAX J1808.4-3658 and 4U 1728-34, and suggest 

£^ | that they are likely strange star candidates. 
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1. Introduction 

The possible existence of a new class of compact stars, which are made entirely 
of deconfined u,d,s quark matter [strange quark matter (SQM)), is one of the 
most intriguing aspects of modern astrophysics. These compact objects are 
called strange stars. They differ from neutron stars, where quarks are confined 
within neutrons, protons, and eventually within other hadrons {e.g. hyperons). 
The investigation of such a possibility is relevant not only for astrophysics, but 
^ ■ for high energy physics too. In fact, the search for a deconfined phase of quark 

matter is one of the main goals in heavy ion physics. Experiments at Brookhaven 
National Lab's Relativistic Heavy Ion Collider (RHIC) and at CERN's Large 
Hadron Collider (LHC), will hopefully clarify this issue in the near future. 

The possibility that strange stars do exist is based on the so called strange 
matter hypothesis, formulated by Witten (1984) (see also Bodmer, 1971). Ac- 
cording to this hypothesis, strange quark matter, in equilibrium with respect 
to the weak interactions, could be the true ground state of strongly interacting 
matter rather than 5e Fe, i.e. the energy per baryon of SQM must fulfil the 
inequality 

E\ E( 5e Fe) 

- < \ } ~ 930 MeV, (1) 

^ J SQM 56 

at the baryon density where the pressure is equal to zero. 

If the strange matter hypothesis is true, then a nucleus with A nucleons, 
could in principle lower its energy by converting to a strangelet (a drop of SQM). 
However, this process requires a very high-order simultaneous weak interactions 
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to convert about a number A of u and d quarks of the nucleus into strange quarks. 
The probability for such a process is extremely low [], and the mean life time for 
an atomic nucleus to decay to a strangelet is much higher than the age of the 
Universe. On the other hand, a step by step production of s quarks, at different 
times, will produce hyperons in the nucleus, i.e. a system (hypernucleus) with a 
higher energy per baryon with respect to the original nucleus. In addition, finite 
size effects (surface and shell effects) place a lower limit (A ~ 10-100) on the 
baryon number of a stable strangelet even if bulk SQM is stable (Farhi & Jaffe, 
1984). Thus, according to the strange matter hypothesis, the ordinary state of 
matter, in which quarks are confined within hadrons, is a metastable state. 

The success of traditional nuclear physics, in explaining an astonishing 
amount of experimental data, provides a clear indication that quarks in a nu- 
cleus are confined within protons and neutrons. Thus, the energy per baryon 
(E/A) u d of u,d quark matter (nonstrange quark matter) must be higher than 
the energy per baryon of nuclei 

(j^j > 930 MeV + A, (2) 

being A ~ 4 MeV a quantity which accounts for the lower energy per baryon of 
a finite chunk (A ~ 250) of nonstrange quark matter with respect to the bulk 
(A — > oo) case (Farhi & Jaffe, 1984). These stability conditions (eq.s (1) and 
(2)) in turn may be used to constrain the parameters entering in models for the 
equation of state (EOS) of SQM. As we show below, the existence of strange 
stars is allowable within the uncertainties inherent in perturbative Quantum 
Chromo-Dynamics (QCD). Thus strange stars may exist in the Universe. 



2. The equation of state for strange quark matter 

From a basic point of view the equation of state for SQM should be calculated 
solving the equations of QCD. As we know, such a fundamental approach is 
presently not doable. Therefore one has to rely on phenomenological models. 
In this work, we discuss two phenomenological models for the EOS of strange 
quark matter. The first one is a well known model related to the MIT bag model 
(Chodos et al. 1974) for hadrons. The second one is a new model developed by 
Dey et al. (1998). 

At very high density SQM behaves as a relativistic gas of weakly interacting 
fermions. This is a consequence of one of the basic features of QCD, namely 
asymptotic freedom. To begin with consider the case of massless quarks, and 
consider gluon exchange interactions to the first order in the QCD structure con- 
stant a c . Under these circumstances the EOS of /3— stable SQM can be written 
in the parametrical form: 

e = Knf + B, P= l -Knf - B, K = ^ (l + ^) he (3) 



1 It is proportional to G 2 ^ , being Gf the Fermi constant, and assuming a number A of simul- 
taneous weak processes. 
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e being the energy density, and P the pressure. Eliminating the baryon number 
density ns one gets: 

P= l -(e-AB) (4) 

Here B is a phenomenological parameter which represents the difference between 
the energy density of "perturbative vacuum" and true QCD vacuum. B is related 
to the "bag constant" which in the MIT bag model for hadrons (Chodos et al. 
1974) gives the confinement of quarks within the hadronic bag. The density of 
zero pressure SQM is just p s = AB / c 2 . This is the value of the surface density of a 
bare strange star. Taking a non-vanishing value for the mass m s of the strange 
quark, the EOS becomes more involved (see e.g. Farhi & Jaffe, 1984) with 
respect to the simple expression (4). However, for m s = 100-300 MeV, equation 
(4) is less than 5% different from the "exact" case for m s ^ 0. In summary, in 
this model for the equation of state for SQM there are three phenomenological 
parameters, namely: B, m s , and a c . It is possible to determine ranges in the 
values of these parameters in which SQM is stable, and nonstrange quark matter 
is not (Farhi & Jaffe, 1984). For example, in the case of non-interacting quarks 
(a c = 0) one has B ~ 57-91 MeV/fm 3 for m s = 0, and B ~ 57-75 MeV/fm 3 for 
m s = 150 MeV. 

The schematic model outlined above becomes less and less trustworthy go- 
ing from very high density region (asymptotic freedom regime) to lower densities, 
where confinement (hadrons formation) takes place. Recently, Dey et al. (1998) 
derived a new EOS for SQM using a "dynamical" density-dependent approach to 
confinement. The EOS by Dey et al. has asymptotic freedom built in, shows con- 
finement at zero baryon density, deconfinement at high density. In this model, 
the quark interaction is described by a colour-Debye-screened inter-quark vector 
potential originating from gluon exchange, and by a density-dependent scalar 
potential which restores chiral symmetry at high density (in the limit of massless 
quarks). The density-dependent scalar potential arises from the density depen- 
dence of the in-medium effective quark masses M q , which, in the model by Dey 
et al. (1998), are taken to depend upon the baryon number density according to 

M q = m q + 310 • sechiu— \ (MeV), (5) 

\ no J 

where no = 0.16 fm -3 is the normal nuclear matter density, q(= u,d,s) is the 
flavor index, and v is a parameter. The effective quark mass M q (ns) goes from 
its constituent masses at zero density, to its current mass m q , as ub goes to 
infinity. Here we consider two different parameterizations of the EOS by Dey et 
al., which correspond to a different choice for the parameter v. The equation of 
state SSI (SS2) corresponds to v = 0.333 (y = 0.286). These two models for the 
EOS give absolutely stable SQM according to the strange matter hypothesis. 

3. Strange star candidates 

To distinguish whether a compact star is a neutron star or a strange star, one 
has to find a clear observational signature. There is a striking qualitative dif- 
ference in the mass-radius (MR) relation of strange stars with respect to that 
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of neutron stars (see Fig. 1). For strange stars with "small" (M << M max ) 
gravitational mass, M is proportional to R 3 . In contrast, neutron stars have 
radii that decrease with increasing mass. This is a consequence of the underly- 
ing interaction between the stellar constituents which makes "low" mass strange 
stars self-bound objects (see e.g. Bombaci 1999) contrary to the case of neutron 
stars which are bound by gravity ^. As we know, there is a minimum mass for a 
neutron star (M m i n ~ 0.1 M Q ). In the case of a strange star, there is essentially 
no minimum mass. As the central density p c — > p s (surface density), a strange 
star (or better a strangelet for very low baryon number) is a self-bound system, 
until the baryon number becomes so low that finite size effects destabilize it. 

3.1. SAX J1808.4-3658 

The transient X-ray burst source SAX J1808. 4-3658 was discovered in September 
1996 by the BeppoSAX satellite. Two bright type-I X-ray bursts were detected, 
each lasting less than 30 seconds. Analysis of the bursts in SAX J1808. 4-3658 in- 
dicates that it is 4 kpc distant and has a peak X-ray luminosity of 6x 10 36 erg/s in 
its bright state, and a X-ray luminosity lower than 10 35 erg/s in quiescence (in't 
Zand 1998). The object is nearly certainly the same as the transient X-ray source 
detected with the Proportional Counter Array (PCA) on board the Rossi X-ray 
Timing Explorer (RXTE) (Marshall, 1998). Coherent pulsations at a period of 
2.49 milliseconds were discovered (Wijnands & van der Klis 1998). The star's 
surface dipolar magnetic moment was derived to be less than 10 26 G cm 3 from de- 
tection of X-ray pulsations at a luminosity of 10 36 erg/s (Wijnands & van der Klis 
1998), consistent with the weak fields expected for type-I X-ray bursters and mil- 
lisecond radio pulsars (MS PSRs) (Bhattacharya & van den Heuvel 1991). The 
binary nature of SAX J1808. 4-3658 was firmly established with the detection of 
a 2 hour orbital period (Chakrabarty & Morgan 1998) as well as with the optical 
identification of the companion star (Roche et al. 1998). SAX J1808. 4-3658 is 
the first pulsar to show both coherent pulsations in its persistent emission and 
X-ray bursts, and by far the fastest-rotating, lowest-field accretion-driven pulsar 
known. It presents direct evidence for the evolutionary link between low-mass 
X-ray binaries (LMXBs) and MS PSRs. SAX J1808.4-3658 is the only known 
LMXB with an MS PSR. 

A mass-radius (MR) relation for the compact star in SAX J1808. 4-3658 
has been recently obtained by Li et al. (1999a) ^ using the following two re- 
quirements, (i) Detection of X-ray pulsations requires that the inner radius Rq 
of the accretion flow should be larger than the stellar radius R. In other words, 
the stellar magnetic field must be strong enough to disrupt the disk flow above 
the stellar surface, (w) The radius Ro must be less than the so-called co-rotation 
radius R c , i.e. the stellar magnetic field must be weak enough that accretion is 
not centrifugally inhibited: 

R <R C = [GMP 2 /(4vr 2 )] 1/3 . (6) 



2 As an idealized example, remember that pure neutron matter is not bound by nuclear forces. 

3 see also Burderi & King (1998), Psaltis & Chakrabarty (1999). 
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Figure 1. Comparison of the mass-radius relation of SAX J1808.4 
-3658 determined from RXTE observations with theoretical models of 
neutron stars and of strange stars. See text for more details. 



Here G is the gravitation constant, M is the mass of the star, and P is the pulse 
period. The inner disk radius Rq is generally evaluated in terms of the Alfven 
radius Rj±, at which the magnetic and material stresses balance (Bhattacharya 
& van den Heuvel 1991): Rq = = ^[B 2 R 6 /M(2GM) 1 / 2 } 2 / 7 , where B and 
M are respectively the surface magnetic field and the mass accretion rate of 
the pulsar, and £ is a parameter of order of unity almost independent of M (Li 
1997, Burderi & King 1998). Since X-ray pulsations in SAX J 1808.4-3658 were 
detected over a wide range of mass accretion rate (say, from M m ; n to M m ax) 5 
the two conditions (i) and (it) give R ^ Ro(M mSuX ) < i?o(-^min) ~ R c . Next, 
we assume that the mass accretion rate M is proportional to the X-ray flux F 
observed with RXTE. This is guaranteed by the fact that the X-ray spectrum of 
SAX J1808. 4-3658 was remarkably stable and there was only slight increase in 
the pulse amplitude when the X-ray luminosity varied by a factor of ~ 100 during 
the 1998 April/May outburst (Gilfanov et al. 1998, Cui et al. 1998, Psaltis & 
Chakrabarty 1999). Therefore, Li et al. (1999a) get the following upper limit of 
the stellar radius: R < (F m i n / F max ) 2 / 7 R c , or 



K < 27.5 ^ C-^fV^-y' 3 km, (7) 

\F ma J \,2A9msJ \M @ J 1 ' 



where F mSuX and F m i n denote the X-ray fluxes measured during X-ray high- and 
low-state, respectively, M is the solar mass. Note that in writing inequality (7) 
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it is assumed that the pulsar's magnetic field is basically dipolar (see Li et al. 
1999a for arguments to support this hypothesis) 

Given the range of X-ray flux at which coherent pulsations were detected, 
inequality (7) defines a limiting curve in the mass-radius plane for SAX J1808.4- 
3658, as plotted in the dashed curve in Fig. 1. The authors of ref. (Li et al. 
1999a) adopted the flux ratio F msx / F m \ n ~ 100 from the observations that during 
the 1998 April/May outburst, the maximum 2-30 keV flux of SAX J1808.4-3658 
at the peak of the outburst was -F max — 3 x 10 -9 ergcm~ 2 s _1 , while the pulse 
signal became barely detectable when the flux dropped below F m j n ~ 2 x 
10" 11 ergcm" 2 s" 1 (Cui et al. 1998, Psaltis & Chakrabarty 1999). The dashed 
line R = R s = 2GM/c 2 represents the Schwartzschild radius - the lower limit 
of the stellar radius to prevent the star collapsing into a black hole. Thus the 
allowed range of the mass and radius of SAX J1808. 4-3658 is the region confined 
by these two dashed curves in Fig. 1. 

In the same figure, we report the theoretical MR relations (solid curves) for 
neutron stars given by some recent realistic models for the EOS of dense matter 
(see Li et al. 1999a for references to the EOS models). Models BBB1 and BBB2 
are relative to "conventional" neutron stars (i.e the core of the star is assumed 
to be composed by an uncharged mixture of neutrons, protons, electrons and 
muons in equilibrium with respect to the weak interaction). The curve labeled 
Hyp depicts the MR relation for a neutron star in which hyperons are considered 
in addition to nucleons as hadronic constituents. The MR curve labeled K~ is 
relative to neutron stars with a Bose-Einstein condensate of negative kaons in 
their cores. It is clearly seen in Fig. 1 that none of the neutron star MR curves 
is consistent with SAX J1808. 4-3658. Including rotational effects will shift the 
MR curves to up-right in Fig. 1 (Datta et al. 1998), and does not help improve 
the consistency between the theoretical neutron star models and observations 
of SAX J1808.4-3658. Therefore SAX J1808.4-3658 is not well described by a 
neutron star model. The curve B90 in Fig. 1 gives the MR relation for strange 
stars described by the schematic EOS (4) with B = 90 MeV/fm 3 . The two 
curves SSI and SS2 give the MR relation for strange stars calculated with the 
EOS by Dey et al. (1998). Figure 1 clearly demonstrates that a strange star 
model is more compatible with SAX J1808. 4-3658 than a neutron star one. 

3.2. 4U 1728-34 

Recently, Li et al. (1999b) investigated possible signatures for the existence of 
strange stars in connection with the newly discovered phenomenon of kilohertz 
quasi-periodic oscillations (kHz QPOs) in the X-ray flux from LMXB (for a re- 
view see van der Klis 2000). Initially, kHz QPO data from various sources were 
interpreted assuming a simple beat-frequency model (see e.g. Kaaret & Ford 
1997). In many cases, two simultaneous kHz QPO peaks ("twin peaks") are ob- 
served. The QPO frequencies vary and are strongly correlated with source flux. 
In the beat -frequency model the highest observed QPO frequency v u is inter- 
preted as the Keplerian orbital frequency vk a t the inner edge of the accretion 



4 see also Psaltis & Chakrabarty (1999) for a study of the influence on the MR rela- 
tion for SAX J1808. 4-3658 of a quadrupole magnetic moment, and of a non-standard disk- 
magnetosphere interaction model. 
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disk. The frequency vi of the lower QPO peak is instead interpreted as the beat 
frequency between uj( and the neutron star spin frequency Uq, which within this 
model is equal to the separation frequency Av = v u — v\ of the two peaks. Thus 
Av is predicted to be constant. Nevertheless, novel observations for different 
kHz QPO sources have challenged this simple beat -frequency model. The most 
striking case is the source 4U 1728-34, where it was found that Av decreases 
significantly, from 349. 3± 1.7 Hz to 278. 7± 11.6 Hz, as the frequency of the lower 
kHz QPO increases (Mendez &: van der Klis 1999). Furthermore, in the spectra 
observed by the RXTE for 4U 1728-34, Ford & van der Klis (1998) found low- 
frequency Lorentian oscillations with frequencies between 10 and 50 Hz. These 
frequencies as well as the break frequency (inbreak) of the power spectrum density 
for the same source were shown to be correlated with u u and v\. 

A different model was recently developed by Osherovich & Titarchuk (1999) 
(see also Titarchuk & Osherovich 1999), who proposed a unified classification of 
kHz QPOs and the related observed low frequency phenomena. In this model, 
kHz QPOs are modeled as Keplerian oscillations under the influence of the 
Coriolis force in a rotating frame of reference (magnetosphere) . The frequency 
V\ of the lower kHz QPO peak is the Keplerian frequency at the outer edge 
of a viscous transition layer between the Keplerian disk and the surface of the 
compact star. The frequency v u is a hybrid frequency related to the rotational 
frequency u m of the star's magnetosphere by: = v\ + (2u m ) 2 . The observed 
low Lorentzian frequency in 4U 1728-34 is suggested to be associated with radial 
oscillations in the viscous transition layer of the disk, whereas the observed break 
frequency is determined by the characteristic diffusion time of the inward motion 
of the matter in the accretion flow (Titarchuk & Osherovich 1999). Predictions 
of this model regarding relations between the QPO frequencies mentioned above 
compare favorably with recent observations for 4U 1728-34, Sco X-l, 4U 1608-52, 
and 4U 1702-429. 

The presence of the break frequency and the correlated Lorentzian frequency 
suggests the introduction of a new scale in the phenomenon. One attractive 
feature of the model by Titarchuk & Osherovich (1999) is the introduction of 
such a scale in the model through the Reynolds number for the accretion flow. 
The best fit for the observed data was obtained by Titarchuk & Osherovich 
(1999) when 

a k = (M/M )(i? o /3i?s) 3/2 (^o/364Hz) = 1.03, (8) 

where M is the stellar mass, R$ is the inner edge of the accretion disk 0, R s is 
the Schwarzschild radius, and vq is the spin frequency of the star. Given the 
364 Hz spin frequency of 4U 1728-34 (Strohmayer et al. 1996), the inner disk 
radius can be derived from the previous equation. Since the innermost radius 
of the disk must be larger than the radius R of the star itself, this leads to a 
mass-dependent upper bound on the stellar radius, 

R<R ~ 8.86 4 /3 (M/M ) 1/3 km, (9) 



5 In the expression for at reported in Titarchuk & Osherovich (1999), one has xq = Ro/R B , where 
Ro is erroneously indicated as the neutron star radius (Titarchuk, private communication). 
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Figure 2. Comparison of the MR relation of 4U 1728-34 determined 
from RXTE observations with theoretical models of neutron stars and 
of strange stars. The range of mass and radius of 4U 1728-34 is allowed 
in the region outlined by the dashed curve R = Rq, the horizontal 
dashed line, and the dashed line R = R s . The solid curves represents 
theoretical MR relations for neutron stars and strange stars. 



which is plotted by dashed curve in Fig. 2. 

A second constraint on the mass and radius of 4U 1728-34 results from the 
requirement that the inner radius Rq of the disk must be larger than the radius 
of the last stable circular orbit i? ms around the star: 

Ro > Rms- (10) 

To make our discussion more transparent, neglect for a moment the rotation of 
the compact star. For a non-rotating star R ms = 3R S , then the second condition 
gives: 

R > 3R S = 8.86 (M/M ) km. (11) 

Therefore, the allowed range of the mass and radius for 4U1728-34 is the region 
in the lower left corner of the MR plane confined by the dashed curve (R = Rq), 
by the horizontal dashed line, and by the Schwartzschild radius (dashed line 
R = R s ). In the same figure, we compare with the theoretical MR relations for 
non-rotating neutron stars and strange stars, for the same models for the EOS 
considered in Fig. 1. It is clear that a strange star model is more compatible with 
4U 1728-34 than a neutron star one. Including the effects of rotation {vq =364 
Hz) in the calculation of the theoretical MR relations and i? ms , does not change 
the previous conclusion (Li et al. 1999b). 
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4. Final remarks 

The main result of the present work (i.e. the likely existence of strange stars) is 
based on the analysis of observational data for the X-ray sources SAX J1808.4- 
3658 and 4U 1728-34. The interpretation of these data is done using standard 
models for the accretion mechanism, which is responsible for the observed phe- 
nomena. The present uncertainties in our knowledge of the accretion mechanism, 
and the disk-magnetosphere interaction, do not allow us to definitely rule out 
the possibility of a neutron star for the two X-ray sources we discussed. For 
example, making a priori the conservative assumption that the compact object 
in SAX J1808. 4-3658 is a neutron star, and using a MR relation similar to our 
eq. (7) Psaltis & Chakrabarty (1999) try to constrain disk-magnetosphere in- 
teraction models or to infer the presence of a quadrupole magnetic moment in 
the compact star. 

SAX J1808.4-3658 and 4U 1728-34 are not the only LMXBs which could 
harbour a strange star. Recent studies have shown that the compact objects 
associated with the X-ray burster 4U 1820-30 (Bombaci 1997), the bursting 
X-ray pulsar GRO J1744-28 (Cheng et al. 1998b) and the X-ray pulsar Her 
X-l (Dey et al. 1998) are likely strange star candidates. For each of these X- 
ray sources (strange star candidates) the conservative assumption of a neutron 
star as the central accretor would require some particular (possibly ad hoc) 
assumption about the nature of the plasma accretion flow and/or the structure 
of the stellar magnetic field. On the other hand, the possibility of a strange star 
gives a simple and unifying picture for all the systems mentioned above. Finally, 
strange stars have also been speculated to model 7-ray bursters (Haensel et al. 
1991, Bombaci & Datta 2000) and soft 7-ray repeaters (Cheng & Dai 1998a). 
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